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Abstract. Our main focus is on a general class of active rotators with mean field 
interactions, that is globally coupled large families of dynamical systems on the unit 
circle with non-trivial stochastic dynamics. The dynamics of each isolated system is 
dipt — —8V'(ipt) dt+ dwt, where V' is a periodic function, w is a Brownian motion and S 
is an intensity parameter. It is well known that the interacting dynamics is accurately de- 
scribed, in the limit of infinitely many interacting components, by a Fokker-Planck PDE 
and the model reduces for S = to a particular case of the Kuramoto synchronization 
model, for which one can show the existence of a stable normally hyperbolic manifold 
of stationary solutions for the corresponding Fokker-Planck equation (we are interested 
in the case in which this manifold is non-trivial, that happens when the interaction is 
sufficiently strong, that is in the synchronized regime of the Kuramoto model). We use 
the robustness of normally hyperbolic structures to infer qualitative and quantitative 
results on the |<5| < Sq cases, with So a suitable threshold: as a matter of fact, we obtain 
an accurate description of the dynamics on the invariant manifold for <5 7^ and we link 
it explicitly to the potential V. This approach allows to have a complete description of 
the phase diagram of the active rotators model, at least for |<5| < So, thus identifying for 
which values of the parameters (notably, noise intensity and/or coupling strength) the 
system exhibits periodic pulse waves or stabilizes at a quiescent resting state. Moreover, 
some of our results are very explicit and this brings a new insight into the combined 
effect of active rotator dynamics, noise and interaction. The links with the literature on 
specific systems, notably neuronal models, are discussed in detail. 
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1. Introduction 

1.1. Coupled excitable systems. There are diverse examples of threshold phenomena 
in natural systems. Dynamics of excitable systems, as exemplified by neuronal membranes 
(to which we restrict for sake of conciseness), constitute one of the common forms of 
threshold behavior. Excitable systems are characterized by their nonlinear response to 
perturbations. In the absence of inputs, they remain at a resting state. This state is 
locally stable in the sense that the system returns rapidly to it after small perturbations. 
However, for inputs beyond a critical range, the response of the system takes on a very 
different form, before regaining the resting state. In the phase portrait of the system, 
subthreshold responses correspond to monotonic returns to the stable equilibrium while 
suprathreshold ones appear as excursions that take the system transiently away from the 
stable equilibrium. Excitability is one of the key neuronal properties at the heart of 
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signal processing and transmission in nervous systems. Motivated by their ubiquity and 
numerous experimental observations attesting to their functional importance, there has 
been a characterization of various forms of excitability in terms of the geometry of the 
phase portrait of dynamical systems [9]. 

Excitable systems are particularly sensitive to noise because such random signals contain 
consecutive sub and suprathreshold segments that occur in an unpredictable manner. The 
interplay between the nonlinearity inherent in the threshold mechanisms and the noise 
induced fluctuations can produce a large variety of dynamics in excitable systems, some of 
which are reviewed in [12]. In this paper, we consider one of these, namely, noise induced 
synchronous coherent oscillations in assemblies of coupled excitable systems. 

Noisy excitable systems display irregular repetitive suprathreshold excursions hence- 
forth referred to as firing. In ensembles of such units receiving independent noise, the 
firings of the units remain independent from one another as long as there are no inter- 
connections between them. Coupling the units with one another introduces correlations 
between their firings. Synchrony is the extreme form of such correlations when the units 
fire almost simultaneously. However, synchronous firings can be irregular. One of the sur- 
prising effects of noise in assemblies of interacting excitable systems is that for some range 
of coupling strength and noise intensity, units fire synchronously and regularly. The wide 
occurrence of these noise induced coherent dynamics and their underlying mechanisms are 
well documented as explained below. Their putative functional role in nervous systems is 
to participate in rhythm generation in the absence of pacemaker units (see for instance 
|10|). Despite the large number of numerical explorations devoted to this phenomenon, it 
has not been analyzed from a mathematical standpoint. The purpose of the present paper 
is to deal with this aspect. 

Two key elements are at play in the occurrence of noise induced regular synchronous 
firing in assemblies of interacting excitable units, one is that interacting excitable units act 
globally like a single excitable system at the population level, the other is that noise driven 
excitable systems undergo coherence resonance [TTJ. How the combination of these two 
phenomena leads to noise induced regularly synchronous firing has been first highlighted 
in an analysis of networks of an elementary neuronal model |18} [T9] , see also [U [20] . 

The important point is the generality of this mechanisms. It neither relies on the refined 
properties of specific classes of excitable systems nor on the types of coupling. In fact, 
noisy assemblies of all common neuronal models, irrespective of the type of excitability, 
and whether coupled diffusively or through excitatory pulses or synapses, readily produce 
noise induced regular synchronous firing. To our knowledge, one of the earliest reports of 
this phenomenon goes back to the explorations of MacGregor and Palasek of randomly 
connected populations of neuromimes incorporating a large array of individual neuronal 
properties [13] . More recent examples include the description of the same phenomenon 
in common neuronal models such as the Hodgkin-Huxley [30], the FitzHugh-Nagumo 
[29], the Morris-Lecar [5], the Hindmarsch-Rose [4] and others implementing detailed 
biophysical properties [10] . In these references, besides differences in the models there are 
also differences in coupling and network architecture: in some the units are diffusively 
coupled, in others they are coupled through excitatory pulses; in some connectivity is all- 
to-all, while others deal with random networks. Our enumeration, which does not intend 
to be exhaustive, illustrates the ease with which assemblies of excitable units generate 
noise-induced synchronous regular activity, irrespective of model and network details. 

The ubiquity of the phenomenon strongly supports investigating its key characteristics 
through the mathematical analysis of a minimal model that captures its essence. The 



TRANSITIONS IN ACTIVE ROTATOR SYSTEMS 



3 



model we consider is a general version of the so-called active rotator (AR) which is repre- 
sentative of the so-called class I excitable systems [9] . 

1.2. Active rotator models. The AR is a variant of the Kuramoto model for excitable- 
oscillatory systems that evolve on a unit circle [T]. Precisely, the AR model can be intro- 
duced via the stochastic equations 

K N 

difjj(t) = -5V (ifjjit)) dt- — J^sin(^(t)-^(t)) dt + adwj{t), (1.1) 

i=\ 

where j = 1, . . . , N, N is a (large) integer, K, a, and 5 are non-negative constants, the 
Wj's are IID standard Brownian motions and V is a smooth function (in the applications 
the case in which V' is a trigonometric polynomial will play an important role, so me may 
as well think of this case). We look at ipj as an element of S := M/2-7rZ, that is ipj is a 
phase, and, of course we have to supply an initial condition for (|l.ip : for example we can 
take {^j(0)}j=i,...,jv to be independent identically distributed random variables. 

This set of equations defines a diffusion on §> N describing the evolution of N noisy 
interacting phases: note that since K > the interaction has a tendency to synchronize 
the ipj's and let us stress from now that such an iV-dimensional diffusion reduces for 5 = 
to a dynamics that is reversible with respect to the Gibbs measure with Hamiltonian given 
by — Yli j cos (4>i ~ V'j) an d inverse temperature a~ 2 . Such a Gibbs measure goes under 
the name of "mean field classical XY model" : we refer to [3] for more details, but we point 
out that for 5 > (of course the case 5 < is absolutely analogous), unless V is a periodic 
function (which we do not assume: consider for example V'^ifj) = 1), the dynamics is not 
reversible. Nevertheless, it is well known that the large N behavior of such a system can 
be described in terms of the Fokker-Planck or McKean-Vlasov PDE (the literature on this 
issue is very vast: see for example the references in [3]): 



dtP S M = ^dU{0)-d e k(0)(J*pf)(0)l +Sd e \p S t (e)V'(9) 



f X2) 



2 

where J(-) := — if sin(-) and & S. To be precise, pf(-) is a probability density and it 
captures the ./V — > oo limit of the empirical (probability) measure J2f=i ^(t)(d^)j where 
5 a is the Dirac delta measure on a. Actually, one can even describe with great accuracy 
(as -/V — > oo) the dynamics of each unit system (in interaction!): it evolves following a non- 
local diffusion equation, called at times non-linear diffusion. The non-locality comes from 
the fact that ipj is subject not only to the force field V, but also to the field corresponding 
to the interaction with all other unit systems, and it all boils down to 

<ty(t) = -5V'(?p(t))dt+{J*pf)('ip{t))dt + (Tdw(t), (1.3) 

with w a standard Brownian motion, and it turns out that the probability distribution of 
tp(t) is precisely pf if ip(0) has distribution p s Q . 

In mathematical terms, the question that we want to tackle is: what is the relation 
between the simple deterministic one dimensional dynamics ip = —V'(i/j) (Isolated Deter- 
ministic System: IDS) and the behavior of the associated N dimensional diffusion, for N 
large? The question is actually twofold. First, given a potential V for the IDS, what is the 
collective dynamic of the N large limit (II. 2D ? Conversely, what are the possible collective 
dynamics of (|1.2p ? In order to be more concrete let us ask the following sharper questions: 
is it possible that 
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• the IDS has only one stable point, for example if V(ip) = ift — acos(tp) for a > 1, 
but the N — > oo system exhibits stable periodic behavior, that is there is a stable 
periodic solution to (|1.2|) ? 

• the IDS has only periodic solutions, but the N — > oo system has stable stationary 
solutions? 

The fact that the answer to these questions is positive is, to a certain extent, known. 
Notably, in their numerical investigations of the dynamics of coupled noisy ARs, Shi- 
nomoto and Kuramoto reported the existence of collective periodic oscillations, the same 
phenomenon we have referred to as noise induced regular synchronous activity [26j . They 
also performed numerical explorations of the transitions to and from this coherent state. 
The key ingredient in such analyses has been to consider the bifurcations of the associated 
Fokker-Planck equation (we anticipate that our results make rigorous some of their pre- 
dictions, see Section [3]). To clarify how noise generates such time-periodic global activity 
in coupled excitable ARs, Kurrer and Schulten approximated the solutions of the nonlin- 
ear Fokker-Planck equations by Gaussian distributions [TT]. Under this assumption, they 
obtained closed ordinary differential equations for the mean and variance of the distribu- 
tion and used the bifurcation diagram of these to investigate the regimes where the model 
generates periodic oscillations. Related work can be found for example in [HHl], where 
finite N analysis has been performed, or in |15l 116} [27], where variants of the model have 
been considered. 

However, from a mathematical viewpoint this phenomenon is only very partially under- 
stood. We are aware of the contributions |21} [22| [23| [28] that are somewhat close in spirit 
to what we are doing: these references deal with periodic behavior in nonlinear Markov 
processes and, more generally, with the effect of the noise on (mean field) interacting dy- 
namical systems. We also deal with nonlinear Markov processes - the evolution equation 
(|1.3I) contains the law of the process itself - even if this aspect is not emphasized in the 
remainder of the paper. In particular, Scheutzow [22J provides examples of mean-field type 
systems in which periodic behavior arises in the N — > oo system, even if it is not present 
in absence of noise. The ingenious model set forth in [22] is however rather particular: for 
example the author plays with some stochastic differential equations of nonlinear Markov 
type that admit also Gaussian solutions and the analysis boils down to studying the be- 
havior of the expectation and covariance of these solutions. This is close to the approach 
taken by Touboul, Hermann and Faugeras [28], who extensively exploit the preservation 
of the Gaussian character that holds for certain nonlinear Markov processes and they do 
so for models that aim at describing neural activity. We stress that in their approach the 
IDS dynamics is linear, while for us the nonlinearity of the IDS is a key feature. Rybko, 
Shlosman and Vladimirov in [21] study a connected network of servers that behaves in a 
periodic fashion in the infinite volume limit, when there are sufficiently many customers 
per server (load per server): in this regime there is an effective synchronization between 
servers and the load per server plays a role which is similar to the parameter K in our 
work, cf. (fTTj) . 

1.3. Informal presentation of approach and results. The purpose of this work is 
to show that for general AR systems one can systematically (at least for some range of 
the parameters) and quantitatively exhibit the relation between the IDS and the infinite 
system. This is done by showing that the (infinite-dimensional!) AR system does behave 
like a one dimensional AR, and the latter can be throughly analyzed. We obtain such a 
drastic reduction of dimension by exploiting the fact that for the 5 = case of (|1.2|) one 
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can perform a rather detailed analysis (due to the fact that it is the grandient flow of a 
free energy functional [3J). In that case and when K > K c := a 2 , stationary solutions of 
(jl.2j) are the constant ^ which is unstable, and a circle M = {q(- — 6q) : 60 G §}, which is 
a manifold of non-constant invariant solutions: these solutions describe the synchronized 
state of the oscillators that have a tendency to be close to 9q. The function q : § — > (0, 00) 
is explicitly known and one can show that M is stable. In fact it has been shown that 
M is stable in the sense that it is a stable normally hyperbolic manifold for the 5 = 
evolution (See Section 12. ip . A deep result in dynamical systems theory guarantees the 
robustness of normal hyperbolicity under suitable perturbations [8] , see also [24] : this 
means that, if 5 > is not too large, there exists an invariant manifold Ms which is stable 
and normally hyperbolic for the evolution (jl.2p . and Ms is a smooth deformation of M. 
In particular, for small enough 5, Ms is still a one dimensional manifold diffeomorphic to 
a circle, and the phase along this manifold plays the role of the natural phase ip G § of 
the IDS ijj = —V'(ip). This makes clearly a direct link between the (one dimensional) IDS 
and the N = 00 system which is an infinite dimensional dynamical system. 

The type of results that we obtain is well exemplified in the most basic of the active 
rotator models, namely the one in which we take V{ip) = ip — acos(ip) (without loss of 
generality: a > 0): note that, for a < 1, the IDS describes just a rotation on the circle, 
while for a > 1 the IDS has a stable point (if) = — arcsin (-)), to which it is driven, unless 
sitting on the unstable stationary point ip = arcsin Q) + 7r. Let us keep in mind that Ms 
is close to M, which is a circle, so that also the dynamics on Ms can be reduced to the 
dynamics of a phase (see Fig. [I]). We are going to show in particular that 

(1) there exists (in fact, we give it explicitly) ao > 1 such that for a G (1,oq) (so the 
IDS has a stable stationary point!) there exist K^,K + > 1, with K_ < K + such 
that for K G (K-, K + ), and 5 > sufficiently small (jl.2p has a stable periodic 
solution - a pulsating wave - which corresponds to the fact that the dynamics on 
Ms is periodic. For K G (1,K—) or K > K + instead the dynamics on the manifold 
Ms has (only) one stable stationary point, so (|1.2|) has a stable stationary solution 
(like the IDS). 

(2) for every a G (0, 1), that is the IDS is rotating, one can find Kq > 1 (sufficiently 
close to 1) such that whenever K G (1,Kq) for 5 sufficiently small the dynamics 
on Ms has (only) one stable stationary point. 



Actually, these examples are just instances of results that we will establish for general 
potentials V. For example we will show that for any V such that V changes sign (so 
the IDS has a stable point), for K large enough, the dynamics on the invariant curve 
stabilizes at an equilibrium for small 5. Or that for any V such that V > (so the IDS 
is rotating), but with nonzero first harmonic coefficient (s), for K close to 1 the dynamics 
on the invariant curve stabilizes at an equilibrium for small 5. 

Finally, regarding the inverse problem, that is the range of possible dynamics, we show 
that given a noise and a coupling strength such that the 5 = system exhibits synchroniza- 
tion, any (phase) dynamics can be produced on Ms, for 5 sufficiently small, by a suitable 
choice of the IDS dynamics (that is, of V) and the relation between these two dynamics 
is explicit. 
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Figure 1. For 8 sufficiently small the solutions of <|1.2[) with an initial condition in a 
L 2 -neighborhood of M (in the figure on the left M is drawn by a dashed line), that is an 
initial condition close to a q^(-) = q(- + i/>), stay close to M for all times. In fact, they 
are attracted by a manifold Ms (solid line, still on the left) that is a small (and smooth) 
deformation of M. For every function q^(-) in M one associates only one function p(-) 
on Ms- While the image on the right stresses the function viewpoint, the one on the left 
stresses the geometric viewpoint: M is a circle and it is hence parametrized just by one 
parameter (the phase ip), but Ms can also be reduced simply to ij). The dynamics on the 
two manifolds is hence reduced to the dynamics of ip, with the substantial difference that 
even if the dynamics for 8 = is trivial in the sense that M is a manifold of stationary 
solutions of (|1.2p with 5 = 0, for 8 > the dynamics on Ms can be non-trivial. As a 
matter of fact, we are going to show that by playing on the choice of V(-) essentially any 
phase dynamics can be observed on Ms, and this for every K and a such that K > a 2 . 



2. Mathematical set-up and main results 
2.1. On the reversible Kuramoto PDE. Let us first sum up a number of results about 

9tP° t (e) = \dfo°M - d e [ P ° t (eKJ*p°tW)] , (2.1) 

where J{6) := — Ksm(8). Note that we have set a = 1: there is of course no loss of 
generality in doing this. We start by introducing the weighted H-\ spaces that are going 
to play an important role in the sequel. 

Given a positive smooth function w : 8 — > (0, oo) we define the Hilbert space H_i^ w as 
the closure of the set of smooth functions S — > R such that J s u = with respect to the 
squared norm || w := J § wU 2 , where U = U w is the primitive of u such that J g wU = 0. 
The alternative way to introduce such a space is in terms of rigged Hilbert spaces can be 
found in [3j. When w(-) = 1 we simply write H_\. Let us remark immediately that 

\H\t 1>Wl = [wi (u W2 - kp^i) < f Wl ul 2 < 

so that the the norms we have introduced are all equivalent. We will also use the affine 
space 

H-! := ue , (2.3) 

provided with the H_i distance. The companion space H\, defined in the analogous way, 
will also appear later on. 



w 2 



-l,W2 



(2.2) 
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Basic features and the stationary solutions of the reversible Kuramoto PDE. The reversible 
Kuramoto PDE has a number of features that we recall here. First of all, the reversible 
Kuramoto PDE has strong regularizing properties [7], so that we can safely talk about 
classical smooth solutions for all positive times, for example whenever the initial condition 
is in L 2 . In particular, (|2.1f) defines an L 2 -semigroup. Actually, the conservative character 
of the dynamics and the fact that we are dealing with probability distributions naturally 
lead to work on the affine space 

L\ := {/ 6 L 2 : // = !}, (2-4) 

with the L 2 distance. One of the main feature of (12. ip . directly inherited from being the 
limit of a reversible stochastic dynamics, is that it is the gradient flow of a free energy 
(which is therefore a Lyapunov functional for the evolution). These properties underlie 
what follows but we do not directly use them, and so we refer to [3], see also [7J for related 
results. 

What plays a direct role in our analysis is the fact that all the stationary solutions of 
(12. ip can be written as 

-exp(2ifrcos(- - ip)) , (2.5) 

Zj 

where ip S S (this accounts for the rotation invariance of (|2.ip ). Z is the normalization 
constant (fixed by the requirement of working with probability densities) and r > is a 
solution of the fixed point problem 

, x As cos(6>) exp(xcos(#)) d6 
r = V(2Kr) with *(x) := Js 1 PV V " . 2.6 

J s exp(j;cos(^)) dO 

^(0) = 0, so that r = is a solution of the fixed point problem and ^- is a stationary 
solution. Moreover ^f(-) is increasing and concave on the positive semi-axis, so that there 
exists at most one positive fixed point r and such a fixed point exists if and only if 
K > K c = 1 (see [3] and references therein). So for K > 1 (that we assume from now on) 
there is a manifold, in fact a curve, of stationary solution, besides the constant solution: 

M := {„(.) = „(. - „) : ■, e S} with ?o (0) := ^^^L. , ,2.7) 

where r = r{K) is the positive fixed point of (|2.6p . We will come back to the manifold 
structure of M, but we point out that the main result in [3] means that M is a stable 
normally hyperbolic manifold ( \24\ p. 494]: we are going to detail this just below): we 
stress that M is actually a manifold of stationary solutions and not only an invariant 
manifold. The key point is that if p® = q G M the linearized evolution operator 

- L q u := i-u" - [u J * q + q J * u}' , (2.8) 

with domain {u £ C 2 (S,M) : f s u = 0} is symmetric in H_ii/ q and its closure, that we 
still call Lg, is a self-adjoint operator operator with compact resolvent, hence the spectrum 
is discrete. Actually the spectrum is in [0, oo): L q q' = and q' generates the whole kernel 
of L q : the spectral gap is therefore positive and it will be denoted Xk (see [3] for a proof 
of all these facts and for an explicit lower bound on A^)- 
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In such a framework it is useful to take advantage of some of the interpolation spaces 
associated to L q . For us the (Hilbert) spaces V q and V q with norms 



\nv q 



'L q v 



(2.11) 

-1,1/? 



V 1 + L q V and IMIv, 2 : = 11(1 + ^)^11-1,1/9: (2-9) 

will play an important role. In [3] it is shown that if v E Lq := {v E L 2 : J s v = 0} 

CiflMb < \\v\\v 9 < c^\\v\\ 2 , (2.10) 

where here (and below) ck denotes a suitable positive constant that depends only on K 
(it will not keep the same value through the text: in particular, in this case it is the same 
for every q £ M). Note that if v E TZ(L g ), 1Z(-) denotes the range of the spectral gap 
guarantees that 

ll*s(n-£ 

At this point it is also worth observing also that, by (|2.2p . there exists ck > such that 
for every ^1,^2 £ § we have 

cjclMI-1,1/^ < ||«||-i,i/^ a < c^lkH-i,!/^ • (2.12) 

Of course, we have the analogous estimates in the case in which l/qy, 2 , or \jq^ x , is replaced 
by 1. 

Stable normally hyperbolic manifolds. We now quickly review the notion of stable normally 
hyperbolic manifold, in the L\ set-up, because it will play a central role in our results. For 
this we need a dynamics: what we have in mind is (|2.20j) but for the moment let us just 
think of an evolution semigroup in L\ that gives rise to {ut}t>o, with uq = u, to which 
we can associate a linear evolution semigroup {$>(u, t)}t>o in Lq, satisfying dt$(u,t)v = 
A(t)<&(u,t)v and $>(u,0)v = v, where A(t) is the operator obtained by linearizing the 
evolution around ut- 

For us a stable normally hyperbolic manifold M C L\ (in reality we are interested only 
in 1-dimensional manifolds, that is curves, but at this stage this does not really play a role) 
of characteristics Ai, A2 (0 < Ai < A2) and C > is a C 1 compact connected manifold 
which is invariant under the dynamics and for every u £ M there exists a projection P°{u) 
on the tangent space of M at u, that is 1Z(P°(u)) =: T U M, which, for v £ Lq, satisfies the 
properties 

(1) for every t > we have 

$(u,t)P°(uo)v = P°{u t )<5>{u,t)v, (2.13) 

(2) we have 

\Mu,t)P°(u )v\\ 2 < Cexp(A 1 t)|| W || 2 , (2.14) 
and, for P s := 1 — P°, we have 

\mu,t)P s {u )v\\ 2 < Cexp(-A 2 0|M| 2 , (2.15) 

for every t > 0; 

(3) there exists a negative continuation of the dynamics {ut}t<o and of the linearized 
semigroup {&(u, t)P°(uo)v}t<o and for any such continuation we have 

\\^(u,t)P o (u )v\\ 2 < Cexp(-Ai*)||t;||2, (2.16) 

for t < 0. 
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As an example - for us a crucial example - let us show that M is a stable normally 
hyperbolic manifold for the L^-semigroup associated to (|2.ip (in Section 12.21 we give some 
details on this semigroup in the general case). As we have seen, M is an invariant manifold: 
it is in fact a set of stationary solutions, so that the dynamics has a (trivial) negative 
continuation, and it is easy to provide an explicit atlas, compatible with the L 2 topology, 
for which M is a C°° manifold and T q M = {aq' : a G R} = K(P°). The projection P° 
we choose is defined by 

t.. .j\ i 

-1,1/9 ' 
W,<l')-l,l/q 

and, since L q q' = for every q G M, we see that Ai can be chosen equal to zero and any 
value C > 1 will do. Moreover if we set v t := Q(q,t)PgV G IZ(Pg) then 



P°(q)v = P°v := , (2.17) 



\vth < c K \\vt\\v q < c K y/l + l/\ 



K 



L q v t 

-1,1/9 

< c K \/l + I/Xk exp(-X K t) Hv^uoll-i^i/g < c' K exp(-X K t)\\v\\ 2 , (2.18) 

where we have used (|2.12|) , then (|2.11|) , then the spectral gap and finally (|2.1U|) . Therefore 
A2 can be chosen equal to Xk, C > c' K , and therefore M is a stable normally hyperbolic 
manifold in L 2 for the reversible Kuramoto evolution, with characteristics 0, Ax and 
C = m&x(c' K , 1). 

For the sequel we observe also that u *— > P£, a map from M to the bounded linear 
operators on Lq, is C°° as it can be easily verified by using for v G Lq the formula 

i, qi , j § J i/ qo 

where, like before, V(9) := JqV, so that V : S — > K is (Holder) continuous and ip \-t 



'IOC 



2.2. The full evolution equation. The type of limit evolution equations we are inter- 
ested in can be cast into the form 

dtP S t (0) = \d 2 p 6 t {e)-d e [p S t {e){J*p 5 t ){9)\ +SG[p 5 t ](6), (2.20) 

where S > and for G we assume 



(1) p 1— > G[p] is a function from L\ to H 



-1. 



(2) there exists 77 > such that G is C 1 (Lf,Lf_i) for every p at L 2 distance at most 
j) of M and the derivative DG is uniformly bounded (in the 77-neighborhood of M 
that we consider). 

Note that p t-t (pJ * p)' is also in C 1 (L^, if-i), in fact even in C°°, so that the evolution 
equation can be cast in the abstract form dp^ = Ap\ + F[pf\ + A complete theory 

of this type of equations can be found in [2JI Ch. 4], in particular for p^ £ L 2 such that 
d L 2(pQ,M) < rj there exists of a unique mild solution in C°([0, T), Lq), for some T > 0. 

Examples include: 

(1) the AR case, that is JT2J, with G[p](0) = d e \p(9)U(9)} and ||E/'|| 00 < 00; 

(2) the case of 

G[p](6) = d e [p(9)J*p(e)] (2.21) 

with J g L°°; 
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(3) the case of 

G[p](9) = deW) [ h(9,9')p(9')], (2.22) 

with h £ L°°, as well as generalizations like dg[p(9) J s h(9,9' ,9")p(9')p(9")] and so 
on. 

In all these examples actually one can prove global well-posedness for arbitrary initial 
condition in L\. But the key point of our analysis is that if the initial condition is suffi- 
ciently close to M, then for 5 smaller than a suitable constant, the solution will stay in a 
neighborhood of M for all times. More precisely, our approach is based on the following 
result, that is essentially contained in [241 Main Theorem, p. 495]. We say essentially 
because the result we need is more explicit for what concerns the various small constants 
that are involved: in Section [5] we detail this issue. 

Theorem 2.1. There exists 5o > such that if 5 £ [0,5o] there exists a stable normally 
hyperbolic manifold Ms in L\ for the perturbed equation (|2.20p . Moreover we can write 

M 5 = + (ps (^) ■ i> E §} , (2.23) 

for a suitable function <ps G C (M,Lq) with the properties that 

• Mq) e n{L q ) ; 

• there exists C > such that sup^(||^5 (q^) H2 + ||$i/,<fo(<2ty,)||2) — C5. 

We are now interested in the dynamics on Mg, which is a curve and, given the mapping 
4>s, the position on the manifold is identified by the phase ipf. A more detailed description 
demands information on nf := (ps{q^s): of course ip® = ipQ and n® = for every t. 

We have the following: 

Theorem 2.2. For 5 G [0,5q] we have that t h-> ipf is C 1 and 

' G W>^)_ 1)1/5 . 

iPf + 5 — *L = 0(5*) , (2.24) 

[Q ,Q )-l,l/q 

with 0{5 2 ) uniform in t. Moreover if we call the unique solution of 

(GfeM) 11/q 

L H n^ = G[qj,] ^ - ' H q^ and = 0, (2.25) 



we have 



supWMVip) - trtyWn, = 0{5 2 ). (2.26) 



A sharper control on the dynamics on Ms can be obtained, under a slightly stronger 
assumption on the perturbation G: it all boils down to go beyond (|2.24p and for this note 
that the left-hand side can be written as R s (ipf) where 



[(ps{q^)J * (psify)]' + $ (G + (ps(qip)} ~ G [q^]) , qL 
RsW := f-r-, • (2-27) 

{q',q)-i,i/ q 

It is clear that Rj, is C 1 , since <p s is C 1 . 



TRANSITIONS IN ACTIVE ROTATOR SYSTEMS 



11 



Theorem 2.3. Under the same assumptions of the previous theorem and assuming in ad- 
dition that DG (recall that G E C 1 (L^; H_\) ) is uniformly continuous in a L 2 -neighborhood 
of Mq ; we have that there exists 5 h-> 1(8), with 1(8) = o(l) as 5 \ ; such that 

sup\R' s (rp)\ < 5 £(5). (2.28) 

V>es 



3. Dynamics on M$: analysis of the active rotators case 

Let us use the results of the previous section to tackle the questions we have raised in 
the introduction for the active rotators case and that, ultimately, boil down to: what is the 
relation between the Isolated Deterministic one dimensional System ip = —V'(ip) (IDS) 
and the behavior of the associated N dimensional diffusion, for N large? So we focus on 
(|2,20p with G[p] = ipV'y and regularity assumptions on V' are going to appear along the 
way. Theorem 12.11 tells us that if || V'Hoo < °o, at least when 5 is small enough, the N — > oo 
limit system - ruled by (12.20P - is described by a dynamics on a one dimensional smooth 
and compact manifold Ms equivalent to a circle and, via Theorem 12.21 and Theorem 12.31 
we have a sharp control on this dynamics. 

In order to be precise on this issue let us speed up time by 1/5 in (12.24p . If we keep 
just the leading terms we are dealing with the dynamics 

i> = -fW), (3-1) 

where / is 

(gmM ... 

W ,Q )-l,l/q 

We say that / 6 C 1 (S,R) - not necessarily the / in (|3.2p - is generic, or hyperbolic, 
if it has a finite number of zeroes on § and all of them are simple, i.e. for all ip for 
which f(ip) = 0, we have f'(il>) 7^ 0- Notice that the set of generic functions is open 
in C 1 (S,M) and dense: if the C 1 distance of / and g is less than (a constant times) e, 
we say that the dynamics generated by / and g are e-close. Note that if e is sufficiently 
small then the two dynamics are topologically equivalent. By this we mean that there 
exists a homeomorphism h : S — > § such that {h(ip(ipo,t)) : t £ R}, where t/)(4>o, •) solves 
ip = —f(tf}) and VKV'OjO) = tpo , coincides with {((j)(h(ipo) , t) : t G R}, where 4>((j>Q,-) 
solves (j) = —g(4>) and 4>((fio,0) = <pQ. Moreover we require that h(-) preserves the time 
orientation, that is for a > sufficiently small and t, \s\ G (0, a] we have that ip(i[)o,t) ^ i/jq 
and h(i/j(ipo,t)) = 4>(h(ipo), s) imply s > 0. 

Theorem 12 . 2 1 and Theorem 12.31 guarantee therefore that for 5 sufficiently small the phase 
dynamics on the Ms manifold speeded up by 

±tf /t = -fU4 /t ) + - § Rs(^ t/s ) , (3.3) 

is (5-close to the dynamics generated by /(•). 

The layout of the remainder of this section is, first, to show that even if we fix K > 1, by 
playing on the choice of V (•), one can generate arbitrary generic phase dynamics on Ms- 
In this part we will make also more explicit the link between V and /. Afterwards, we will 
work out in detail a few particular cases and expose some a priori surprising behaviors, 
notably that IDS with periodic behavior (active state) may lead to a iV — > oo dynamics 
that settles down to a fixed point (quiescent state) or that IDS without periodic behavior 
may give origin to periodic iV — > oo behaviors. 



12 



G. GIACOMIN, K. PAKDAMAN, X. PELLEGRIN, C. POQUET 



3.1. Noise and interaction induce arbitrary generic dynamics. It is practical and 
sufficient to work with V'(-) that is a trigonometric polynomial, that is 

n 

V'{9) = a + ]T ( aj cos(j9) + bj sin(j0)) . (3.4) 

3=1 



Theorem 3.1. For any generic dynamics on the circle ipt = —f(ipt) with f £ C 1 (S;IR) 
and for any value of K > 1 there exists a trigonometric polynomial V'(-) (see Remark \3.2\ 
for an explicit expression) such that for 5 small enough, the phase dynamics on Ms (j3.3[) 
is 5-close to ip = — flip). 

Proof. Let / be a generic function in C 1 . By the Stone- Weierstrass Theorem, for every 
e > there exists a trigonometric polynomial P(-) such that \\f — P\\oo < e - If c o is such 
that fn(P — Co) = then, since J Q n f = 0, \cq\ < £■ Thus if we define the trigonometric 
polynomial Q(ip) := /(0) + $(P(9) - c ) d9 we have 

HQ- /lid = \\Q ~ /Hoc + \\P -co -/'Hoc < (27r + l)||P-c -/ / ||oo < (4vr + 2)e, (3.5) 
so it suffices to consider functions / which are trigonometric polynomials: 

n 

f(8) = Ao + Y, ( A k cos(k9) + B k shx{k9)) . (3.6) 
fe=i 

Now we observe that if V'(-) is of the form (I3.4h then a straightforward calculation gives 



(?', Q')-i,i/ q 



k=l 

where 

r-2-K 



I 

a ° + 72 — 7 ^2 ( IkCLk cos (^) + J kh am(kip)) , (3.7) 



I k = I k {2Kr{K)) := — [ * cos{k9)e 2Kr{K)cos{ - e) d9 . (3.8) 
2vr Jo 

Therefore by making the choice ao := Aq and for k = 1, 2, . . . , n 

I 2 - I I 2 - 1 

a k := -2—- A k and b k := ~^-^B k , (3.9) 

we obtain the function V'(-) we were after. □ 



Remark 3.2. The link between / and V can be made more explicit. In fact from ([3? 
and (|3.9p and the fact that the Fourier series of qo is 



»<*> " 2^2K7/ K "°- m = h 4 | W) ""<*> • (3 ' 10) 
one directly extracts that 

/ = ao + j^^r- 1 (9 ° * ^ " ao) = a ° + D{K)q ° * (y ' " ao) ' (3 - U) 
where we have set 

. . I%(2Kr(K)) 



TRANSITIONS IN ACTIVE ROTATOR SYSTEMS 



13 



and (|3.11|) can be applied also in the case in which / is not a trigonometric polynomial. It 
tells us that, for S small, the effective force that drives the N — > oo system is, in a sense, 
obtained by smearing V' via the probability kernel qo. To be precise, V — ao is smeared 
and multiplied by D(K), while the th order Fourier coefficient is left unchanged. This is 
telling us that the effect of noise and interaction, to leading order, boil down to the size of 
D(K) and to the smearing effect of the probability kernel qo(-) (that depends on if too!). 

While (|3.1ip is quite explicit, it is not always straightforward to read off it the qualitative 
properties of /. We start by analyzing the case of if very large and the case of if close to 
one, before moving to treating in detail some particular cases. 

The if — > oo limit. It is straightforward to see that the probability density qo(-) con- 
verges to the Dirac delta measure at the origin. Moreover lmiR-_ il00 D(K) = 1, since 
lmiR-^oo r(if) = 1 and lim x _ ) . 0O Iq (x) = oo. Therefore / and V get closer and closer as 
if becomes large. More precisely one has that for every s £ N and every trigonometric 
polynomial V'(-) there is C such that 

11/ - Vile < -j=. (3.13) 



The proof can be obtained for example by using f|3.9|) that, with (|3.12p . tells us that Aj/aj, 
as well as Bj/bj, that is the ratio of the (non- vanishing) sine and cosine Fourier coefficients 
of / and V , is 

I (2ifr(if)) 

so that by using (Ij(x)/Io(x)) — 1 x ^°° |- (J = 1,2,...) and linifc;_ s . 0O D(K) = 1 we 
readily obtain that the j th -Fourier coefficients of /(•) are, to leading order, j 2 aj/(4if ) and 
j 2 bj /(4K). Since we are just dealing with trigonometric polynomials and the estimate of 
the 1? norm of arbitrary derivatives of f — V, via Parseval formula, is straightforward, 
we get to (|3.13|) . This means in particular that given a potential V such that V' has sign 
changes (so that the IDS has stable points), for any K large enough, the N — > oo system 
has stable stationary solutions, for 5 small enough. We will encounter this phenomenon 
in the particular cases that we treat below. 

The if \ 1 limit. This time we use r (K) ^/2{K - 1) and we derive, first of all, 

that D(K) ~ (4(if — 1)) _1 , since Iq(x) — 1 x 2 /4. Once again we analyze the Fourier 
coefficients of /, via (|3.14p . and we use for j = 1, 2, . . . 

A i = B i Ky (g-l)-i+WP 
a 3 bj 2 2 -0'/ 2 )j! ' 1 ] 

Notably, the first Fourier coefficients of / are enhanced with respect to the corresponding 
coefficients of V by a factor that diverges like (K — 1)~ 1 / 2 . The second Fourier coefficients 
of / are (asymptotically) just proportional to the ones of V' , while higher coefficients in 
the if \ 1 limit are depressed passing from IDS to N — > oo behavior (recall that the 
th -order coefficient is unchanged) . A quantitative estimate in the spirit of (|3.13|) is easily 
established from these estimates. 



so that for j = 1, 2, 
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What we retain from this K \ 1 analysis is that if the first Fourier coefficients are 
present, that is |oi| + \b%\ > 0, then for K sufficiently close to one f(ip) = has two 
solutions and the dynamics will eventually settle to a fixed point (quiescent state). If 
instead |oi| + = 0, then it depends on the relative size of ao and 02 or 62 whether the 
system is in an activated or quiescent regime. But if also 02 = 62 = (and ao 7^ 0) then 
for K sufficiently close to one we have that f(ip) is close to ao and therefore f(4>) 7^ for 
all V>> so that the dynamics is periodic. Again, we will discuss in more detail these issues 
below, in specific examples. 

Remark 3.3. The analysis for K large and close to one is helpful to get an idea on the 
relation between / and V , but the reader should keep in mind that the 5-closeness of the 
dynamics holds for fixed K, that is for 5 < 5q(K). Quantitative estimates on how 5q(K) 
behaves for extreme values of K is an interesting issue that we do not approach here. 

3.2. Active rotators with V(9) = 6 — acos(9). Without loss of generality we assume 
a > 0. Let us start the analysis by making a remark on the a = case: the potential 
becomes just a straight line, and (|2.20p reads 



In this case pf(9 — St) solves (|2.ip . thus Ms = M and the dynamics on Ms is a rotation 
for all 5. 

If a > we exploit the analysis we have developed for Theorem 13.11 that tells us that 
the N — > 00 phase dynamics is lead by the effective force 



Therefore if a < a c (K), then the dynamic on Ms is periodic for 5 small enough ( depending 
on K ) and if a > a c (K), there are two fixed points. From this observation and the graph 
of a c (-) (see Figure [2]) we draw the following conclusions (see also Figured]): 

• Set a c := max^ a c {K){> 1). If a > a c then for every K we have that f{6) = has 
two solutions, so that the phase dynamics has two stationary hyperbolic point: one 
is stable and the other is unstable. In this case the dynamics of the IDS resembles 
to the phase dynamics of the N — > 00 system. 

• If a & (1,6c) then a c (K) = a has two solutions AT_(a) < A^ + (a) and for K £ 
(K_(a), K + (a)) we have a < a c (K), that is f(9) < for every 9, and the motion 
is periodic: in this case the dynamics of the IDS, that has two fixed points, differs 
from the N — > 00 phase dynamics. For K > K + (a) and for K < K-(a) instead 
the phase dynamics is driven to a (unique) stable fixed point (unless it starts from 
the unstable fixed point). 

• If a < 1 instead a c (K) = a has only one solution K{a) and the periodic behavior 
sets up for K > K(a), otherwise {K < K{a)) the system eventually settles on a 
fixed point: this second case is another instance in which the dynamics of the IDS 
and the N — > 00 system differ. 

When the phase dynamics is periodic we can explicitly integrate the evolution equation 
(|3.ip and compute the first order approximation the period Ts(a, K) of the dynamics on 



d tP s t(9) = ^d 2 eP 5 t (e)-d e [p 5 t (e)(J*p 5 t )(9)-5p s t (e) 



(3.17) 




(3.18) 



M s : 



T s (a,K) 



5 



+ O(l) , where r(a, K) : 



VI - {a/a c {K)Y ' 



(3.19) 
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Figure 2. The graph of a c (-). For K -* oo we have a c (K) = 1 + 1/(8K) + 0(K~ 2 ), 
while for K \ 1 we have a c (K) ~ N /32(iif - 1). 



Actually, it is possible to replace in this formula 0(1) with 0(5): in fact it is possible to 
show by induction that the phase speed on Mg admits an expansion in (integer) powers of 
5 to any order (but with coefficients less explicit than the first order one), and it is easy 
to see that ipg is an odd function of 5. We have tested numerically this approximation and 
we report the result in Table HJ 




a c (K) a 

Figure 3. A sketch of the phase behavior for V(6) = 6 — acos(6): for a > K there 
are two fixed points, one attractive and one repulsive, while for a < a c (K) the force is 
bounded away from zero and the motion is periodic. 
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Table 1. We have simulated (|1.2[) with V(9) = 
9 — acos(8) for a — 1.1 and K = 2. In this 
case our estimates ensure the existence of peri- 
odic solutions for 8 sufficiently small and the pe- 
riod given in (|3.19p (in fact, r := r(l.l, 2) = 
18.0779...). The simulation, that has been 
performed via Fourier decomposition (50 modes 
kept), gives c = 0.333 . . ., for the constant c such 
that (£T 4 (1.1,2)/t(1.1,2)) - 1 ~ c8 2 . 



5 


ZKl.1,2) 


t(1.1,2)/<5 


0.005 


3615.59 


3615.62 


0.010 


1807.79 


1807.85 


0.020 


903.89 


904.01 


0.040 


451.94 


452.19 


0.080 


225.97 


226.45 


0.160 


112.98 


113.96 


0.320 


56.49 


58.51 


0.640 


28.24 


33.02 



3.3. Active rotators with V(9) = 6 — acos(jO)/ j , j = 2,3, . . .. In this case the N — > oo 
phase dynamics is lead by 



1 + a J 2°J ? 1 sin C^)) > 



(3.20) 



and the behavior differs substantially from the j = 1 case (and the j = 2 case is different 
from the j > 3 case). In this case the crucial function is 



a c ,j(K) :-- 



i o 



(3.21) 



Note that a Cj i = a c . The criterion to have periodic behavior is, like for the j = 1 case, 
a < a c j(K), while a > a c j(K) leads to two fixed points. Figure [Hand its caption describes 
the (relatively surprising) phenomenology of the j = 2 and j = 3 cases (the case j > 3 is 
qualitatively the same as the case j = 3). 





Figure 4. For V(0) = 9 — acos(j9)/j, j > 2, the TV — ► oo dynamics is always periodic 
for a < 1 (and 8 sufficiently small), unlike the j = 1 case (recall that a < a c j(K) 
corresponds to periodic motion, while a > a Ci j(K) corresponds to two fixed points: see 
the text). Moreover, for j = 2 and a > 4 the dynamics has just two fixed points, but for 
j > 3 for arbitrarily large values of a one can observe periodic motion if K is sufficiently 
close to 1 (and, of course, 8 sufficiently small). 
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Remark 3.4. Theorem 13.11 already tells us that one can produce arbitrary dynamics, 
so a very large variety of phenomena is observed. Here is a case that can be of some 
interest since it shows that playing on only one parameter one can produce three different 
dynamics (and the reader will directly infer how to induce arbitrarily many): if V{6) = 
9 — a(cos(#) + cos(2#)) the N — > oo phase dynamics is lead by 

f($) = 1 + a-fi- (h sin(V) + 2/ 2 sin(2^)) , (3.22) 

and in this case there can be two transitions as a varies. For example for K = 2 we have 
periodic behavior for a < 0.600 . . ., two fixed points if a € (0.600 . . . , 2.107 . . .) and four 
fixed points (of course two stable and two unstable ones) if a > 2.107 .... 

4. Perturbation arguments 

In this section we assume that 5 G (0, 5q] (cf. Theorem 12. ip and that we are on the 
invariant manifold Ms of (I2.20p . that is pf G Ms for every t. The result [241 Main Theorem, 
p. 495] actually contains also some estimates on the regularity of the semigroup on Ms and 
notably that 1 1— > p\ belongs to C°(M;Hi) and that it is (strongly) differentiable as a map 
from M. to H—\. One directly sees that \\u — v\\h-l = \\u" — t/'||-i, so that the right-hand 
side in ([2T20]) is C°(K; and, in turn, t ■->• p s t is C7 1 (R, H-i). 

Since we are working in a neighborhood of M it is useful to introduce from now a 
parametrization of this region that will be particularly useful in the next section, but that 
we are going to use from now. The following facts are proven in Lemma 15. II for every u 
in a sufficiently small H-\ neighborhood of M there exists a unique q = v(u) € M such 
that 

(u-q,q')_ 1A/q = 0. (4.1) 

Furthemore v G C 1 (-fT_i, i/_i) with differential 

Dv(u) = P° (u) . (4.2) 

Theorem 12.11 is telling us in particular that 

v(q + fo{q)) = q- (4.3) 

For the arguments that follow it is practical to use the notation introduced right after 
Theorem 12.11 and write 

pf = q^+nt (4.4) 



where 



q^ t = v (p s t) and n t ■= 4>& \q^i) ■ 



Proof of Theorem {KM Since the evolution on Ms is C^M, H-\), then t i— > q^s is 
C 1 (M, H-i) too. This implies that, with f\ and f2 respectively sine and cosine, ip h-> 
f s q 1 p(9)fi(6) =: aj(t) is C 1 . Since fi(ipt) = °« (t) I V a i (*) + °2 (*) > we see * na ^ ^ ^ is 
C 1 . The fact that and ^I 5 are C 1 directly implies that i i— >• is C 1 (R; i?-i) (actually, 
since ^ is C 1 we have even n s G C 1 (M; Lq)). 
Notice furthermore that 

-Mi = po %M- ^ 

This follows by taking the time derivative of both sides of the equality q^s = v(pf) and by 
using (|4T2"j) . 
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Using (|2.2U|) and the fact that q^, is a stationary solution of (|2.1|) . we rewrite (|4.5p as 



SJ 



-d e 



4(J*n S t ) 



+ SG 



Recall that ||n t 5 || 2 < CS (cf. Theorem EU) : by 



nU * nf 



-i < \\Jh\\n S t \\ 2 2 < C 2 \\J\\ 2 S 2 , 



and by the hypothesis on G that implies that 



G 



from (14.61) we see that 



m + 



G 



5G 



< c G C5. 



< cS 2 , 



(4.6) 

(4.7) 
(4.8) 
(4.9) 



with c independent of t and of ipQ. To obtain (|2.24p just take the H-i >q s scalar product 
of q' s and the expression inside the norm in the left-hand side of (|4.9p . 



For (p^6D rewrite (j2T20|) as 



6 J 



HI 



d t n\ 



+ 5G 



n, 



Note that for the second term on the left hand side we have 



where we have use 

and the bound on the derivative of <j>g given in Theorem 12. 11 

Now plug (pOlj) into (jUniD and use (gTfJ), (gSJ) and KTT\i to obtain 



sup 

t^ s 



( 



G 



G 



(4.10) 

(4.11) 
(4.12) 



\ 



V 



(gW)-l,l/g 



HI 



J 



Q{S 2 



(4.13) 



Since can be chosen arbitrarily on S, we can replace ipf with ip and take the supremum 
over ip (and, by (|2,2p . we can freely switch between H_\ and H_x,\/q^, norms). Therefore 
(recall ([2T25D ) 



sup 



0(<5 2 ) . 



(4.14) 



There result we are after, that is (|2.26p . follows from the equivalence of H\ and VI 2 (recall 



(|2.9p ) norms, which is proven in Appendix [Aj □ 

Proof of Theorem \2.3[ It is of course sufficient to estimate the numerator in the right-hand 
side of ()2.27p . It is the sum of two terms: the first one can be rewritten as 

and, by derivating and using the two L 2 -estimates on (f>s{-) and D<f> s (-) in Theorem 12. 1\ it 
is straightforward to see that there exists c > such that for S £ [0, So] 



(4.15) 



sup|T{(VOI < cS 2 



(4.16) 
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Let us turn to the second term, that is 

W) = jf (1 - 2n j\f q e) ) [ (g[**W) + tfwe'))} - GfaW]) de'do. (4.17) 

For this we write 

H[y] = G [y + ^(y)]-G[y]. (4.18) 

We have 

DH[y] = DG[y + <j>\y)\ - DG[y] + DG[y + <t>\y)]D<j> s (y) (4.19) 

and thus, using the estimates of theorem l2.1l and the fact that DG is uniformly continuous 
on a neighborhood of M, we get that 

sup|T^)| < 1(5). (4.20) 

V>es 

with 1(5) = o(5) when 5 — > 0. 

□ 

5. On the persistence of normally hyperbolic MANIFOLDS 

In this section we prove theorem 12.11 The proof in a more general case can be found in 
[24] but we pay more attention on the relation between the various small parameters that 
enter the proof. We first give a lemma which defines a parametrisation in a neighbourhood 
of M using the scalar structure given by the operators L q . The proof of this lemma is in 
[211 p. 501]. 

Lemma 5.1. There exists a a > such that for all p in the neighborhood 

N a := U qeM B L 2(q,a), (5.1) 

of M there is one and only one q = v(p) G M such that (p — q, <3 ,/ )-l,l/9 = ^' Furthermore 
the mapping p i— > v(p) is in C 00 (L\,L^), and 

Dv(p) = P° v[p) . (5.2) 

Moreover, the analogous statement holds if N a is replaced by U g gM-Bf/_i (q, o~) and this 
time p ^ v(p) is in C°°(H-\,H-i). 

For the proof we look for conditions on 5 in order to get a manifold, which is invariant 
for for (|2.20p . at distance e from M: the condition in the end is going to be that 5 needs 
to be smaller than a suitable constant times e (and e sufficiently small too), so that the 
invariant manifold is in a neighborhood of order 5 of M. To simplify notations, we will 
write F[u] = dg(uJ * u), and ()2.20p becomes: 

dm = \d 2 e pt - F]pt] + SG[ Pt ] . (5.3) 

We will consider solutions with initial condition p satisfying \\p — ^ (po ) 1 1 2 ^ £• We need 
asumptions on e and 5 such that the solution stays in N a for a sufficiently long time. If q 
is in M, w t ■= Pt — q satisfies 

wt = e - tL *w + I e - (t - s)L « (F[w s ] + 5G[q + w s }) ds , (5.4) 
J o 

and we get 

1Mb ^ IKH2+ ri|e- (i ~ s)i9 ll£(H_ 1 ,L 2 )(l|i ? Mlk_ 1 +<J||G'[g + U ; s ]||^_ 1 )d S . (5.5) 
J 
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Define 

to = sup{i ^ : 1 1 i/7 s 1 1 2 ^ c for every s ^ t} . (5-6) 

Because of the continuity of wt, to > if we suppose e < a. If t ^ to, using the spectral 
properties of L q and the regularity of F and G, we get the bounds 

l|e- (t - s)L * \\c(h. u l 2 ) <C L (l + (t- S )-V2) , (5.7) 

\\G[q + w s ]\\ H ^ ^ C g (1 + |K|| 2 ), (5.8) 

and 

||*K]||ir-i < C F \\w s f 2 , (5.9) 

and thus for all t\ < to 



\w tl h ^ (s+C G C L (t 1 + 2^)S) + C L (C F a+C G 6) (l + — L 



\w s \\ 2 ds. (5.10) 

We need the following lemma, that is a version of the Gronwall-Henry inequality 

Lemma 5.2. Let t ^ y t be a non-negative and continuous function on [0,T) satisfying 
for all t e [0, T) 



yt < % + m y o ( i + -7= J y s ds . (5.ii) 

T/ien for all t € [0, T) 

Vt ^ 2r] e at , (5.12) 
with a = 2r]i + (T (|)) 2 where T(r) = J °° x r_1 e _:r dx. 

Proof of lemma \5.2\ We consider the time 

t* = sup{t ^0,y s ^ 2r] e as for all s t} . (5.13) 
We have to show that t* = T. But if t* < T, then 

y t * ^ m ( 1 + 2r/! / ( 1 + -yp=) «s j 



< % (l + ^V'* - 1] + Q j e-j < 2noef* , (5.14) 

which contradicts t* < T since y. is continuous. □ 
Using Lemma 15.21 and (|5.1U|) we get : 

IKIh < C(h){5 + e), (5.15) 

where 

Cifx) = max(l, C G C L (h + 2 ^ti')) e ( 2 '^>< 5 )+W-,<5) 2 >i , (5.16) 

r,(a,6) = C L (C F a + C G 5). (5.17) 

For T > 0, if we choose e and <5 such that C(2T)(e + 5) ^ cr, then lies in iV CT for 
t £ [0; 2T]. Take now T such that 

C ps e-^ < I, (5.18) 
lb 

e AlT/2 ^ 4C L , (5.19) 
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where we recall that Ai is the spectral gap of L q and we set 

C P s = max||P 9 s || £(L g >Lg) , (5.20) 

and Pq is a compact notation for P s (q) (defined just below ()2.14[) ) and it is the orthogonal 
projection of the range of L q (the scalar product is the one of H-\ t \/q)- Define also 

Ci = C(2T), C 2 = e AlT / 2 and e = . (5.21) 

For now we will take max{e, 5} ^ So, so that pt G N a for t < 2T. We will use the following 
notations: 

Pi := p(t, Pi0 ), (5.22) 

is the solution of f)5.3j) and 

Vi := Vi(t,pio) := u(jPt), (5.23) 
is given by Lemma 15. II Moreover we set 

"i =Pi~Vi, Ap := p\—p 2 , Av := vi - v 2 , An := n x - n 2 . (5-24) 

In the following lemma we compare the quantites we have just introduced with the initial 
conditions. It corresponds to Lemma 74.7 (page 507) in |24j . We remark that is in this 
lemma e and 5 play the same role and we stress that these are just preliminary estimates: 
some of them are going to be refined later on. 

Lemma 5.3. For all a > 0, there exist Co = Cq{T) and e\ ^ eq such that if £ ^ e\ and 
5 ^ Ex we have the following properties: 

(1) if \\Po ~ v \U < e then for all t G [0, 2T] 

max (\\p(t,po) - v h , \\v{t,p ) - v \\ 2 , ^||n(t,p )||2^ ^ C (e + S) ; (5.25) 

(2) if\\pio - Vioh < £ and ||Av(0)|| 2 < ae, then for all t G [0, 2T] 

max(||Ap(t)|| 2 , ||Aw(*)|| 2 , ||An(i)|| 2 ) < C 2 ||Ap(0)|| 2 , (5.26) 
with C 2 given in (15.211) ; 

(3) i/|bi0-u i0 || 2 < e and ||Ap(0)|| 2 2||Av(0)|| 2 , then for all t G [0,2T] 

~||A W (0)|| 2 < ||Aw(t)|| 2 < |||A«(0)|| 2 . (5.27) 

Proof o f Lemma 1 5. 3\ For what concerns part (1) note that the first of the three inequalities 
in (]5.25p is given above (see (|5.15p with to = 2T). The other inequalities come from the 
fact that the mapping q \— > v(q) of Lemma 15.11 is Lipschitz, taking, if necessary, a bigger 
value for Co- 

For part (2)notice that, since v 2 q G M, we can write the evolution in mild form around 
v 2 q, that is 

Ap(t) = e~ tLv 2oAp(0) 

+ f e-(*- s ^o (F\px(a) - v 20 ] - F[p 2 {s) - v 20 ] + 5(G[ Pl (a)] - G[p 2 (s)])) ds , (5.28) 
J o 

and thus 

||Ap(t)|| 2 < C L ||Ap(0)|| 2 + C L (c F (ae + C (e + <5) + C G 5)^l + ^=L=^) ||Ap(s)|| 2 ds . 

(5.29) 
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Here we used the preceding point, (|5.7[) and the bounds 

\\G[ Pl (s)]-G\p 2 (s)}U < C G || Pl ( S )-p 2 ( S )|| 2 , (5.30) 

||FbiW]--F[P2W]||-i *S C F ( a e + C (e + 5))\\p 1 (s)-p 2 (s)\\ 2 , (5.31) 
(|5.3ip is obtained by applying the mean value inequality to F and DF and using the fact 
that DF(0) = 0: the constants Cq and Cf have a larger value than in (j5.8j) and (|5.9|) . 
Applying Lemma 15.21 to (|5.29p . we obtain 

||Ap(t)|| 2 < 2C L e( 2 " 1 ( e ' <5 ) +4 ^ 1 ( e ' 5 ) 2 ) 2T ||Ap(0)|| 2 (5.32) 

with 

r/i(e,5) = C L (c F {ae + 2C {£ + 5)) + C G 5) , (5.33) 
Choose £\ ^ Eq such that (it is possible because of (|5.19p ) 

2CLe( 2m(£l ' £l)+47rm(El ' El)2 ) 2T ^ e AlT/2 . (5.34) 

The two other points come directly from the Lipschitz property of the mapping q \— > v(q) 
taking, if necessary, a smaller value for e\ . 

For part (3) we prove first that for all r > 0, there exists £2(f) such that for all e sj 62(f) 
and 5 ^ £ 2 (r) we have for all i E [0, 2T] 



1 ||Au(t)|| 2 3 

2 < k» * 2 (5 - 35) 



if ||Av(0)|| 2 r. In fact, in this case, using Lemma 15. II : 
|A«(i)|| 2 -||Ai;(0)||3 



|A<0)||, 



^ I ||A«(t)|| 2 -||A^0)|| 2 I 
r 

||A»(.)-A.,(0)|h 
r 

< \\vi(t) ~ Vl(0)h + \\v2(t) ~ V2(0)h 

2C (5 + e) 
r 

We can choose £ 2 (r) = min(ei, t/SCq). Now it is sufficient to prove that, for ||Au(0)|| 2 ^ tq 
with a certain ro, 

||Au(i)-A«(0)|| 2 < ^||Aw(0)|| 2 (5.37) 

for all t E [0,2T]. Suppose that ||Au(0)|| 2 ^ r with r ^ a. We use the following decom- 
position 

Av(t) - Av(0) = Av(t) - P° 2 Ap(t) - Au(0) + P° Ap(0) (5.38) 
+ (P° - P° )Ap(t) + P° (Ap(t) - Ap(0)) . 

From Lemma I5TT1 . part (2) and the hypothesis ||Ap(0)|| 2 ^ 2||Av(0)|| 2 , we get 

\\Av(t) - P° 2 Ap(t)\\ 2 = \\v(pi) - v( P2 ) - Dv V2 (pi - p 2 )\\2 

< 2C t ,(r + 2C7o(5 + e))||Au(0)||2 (5.39) 
||Au(0) - P° 20 Ap(0)|| 2 = \\v(p w ) - v(p2o) - Dv V20 (pi - P2o)h 

< 4C v r\\Av(0)\\ 2 (5.40) 
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where C v is the maximum of the second derivate of q i— > v(q) in N a . Since P° is C 1 on 
M, there exits Lpo such that (By using part (1) and the hypothesis and defining Lpo the 
maximum of the nor me of DP on M, we get 

\\(P v ° 2 -P v °JAp(t)\\ 2 ^L P oC (5 + e)\\Av(0)h (5.41) 
For the last term, we write 

Ap(t) - Ap(0) = (e~ tL ™ - J)Ap(O) + f e-W*** (F\p ls - v 20 ] - F[p 2s - v w ] 

Jo 

+ G\pu] ~ G[ P2s ])ds . (5.42) 
Notice that P° 2o (e~ tL ^o - J)Ap(O) = 0. From ([521), flEHJ) and (px30l) it comes 

||P° (Ap(t) - Ap(0))|| 2 < 4C P oC7 L C 2 (c F (r + 2C7 (<5 + e)) + C G 5) (T + v/2f)||A«(0)|| 2 

(5-43) 

where Cpo is the maximum of the norms \\P° \\c(l 2 ,l 2 ) f° r 9 ^ I 11 conclusion there 
exists a constant C3 such that for all t € [0, 2T] 

||Au(i)-Au(0)|| 2 < C 3 (r + e + 5) ||Av(0)|| 2 . (5.44) 

To end the proof choose r = ro := min (a, ^) and reduce if nececessary the value of e± 
to have E\ ^ min(e 2 (ro), ro). 

□ 

We now move to the main body of the proof which is based on introducing a family of 
transformations of the manifold M by using the full dynamics and we aim at identifying 
the transformation that maps M to the manifold that is stable for the full dynamics 
and this is achieved by applying the Banach fixed point Theorem in a relevant space of 
functions. 

Define the set C(M,Lq) of continuous functions from M to L\ provided with the norm 



» = sup{||/(u)|| L 2 , v€ M} (5.45) 

and consider the subset F(e,l) of C(M,Lq) of functions / satisfying : 

(1) ll/lloo<e 

(2) / is Lipschitz on M with Lipschitz constant / ^ 1 

(3) (/(<?), = for all gin M 

Notice that T{e, I) is a complete subset of C(M, Lg). We will now define a set of mappings 
{X t } te r T 2T] '■ F{ £ i 1) ^ C(M,Ll) and show that 

(1) for all r G [T,2T] 

X T (f(e,l)) c?(e,-^\ (5.46) 

(recall that C 2 = e XlT/2 and thus ^ ^ 1) 

(2) Xt is a contraction on J-(e, 1): 

\\X T {h) - *t(/ 2 )IU ^ \\\h - /alloc (5-47) 
for all /i,/ 2 € T(e,l). 
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Notice that the third point of (|5.3|) and an argument of connexion ( see |24j page 513 ) 
show that for all / G T(e, 1) and t G [0, 2T], the mapping q h-> gtj(q) '■= v(t,f(q)) is a 
bijection of M. So we can define the mappings 



^ n(r,(t d + /jo (5 ' 48) 



It is easy to see that for all r G [T, 2T] and / G .F(£, 1), X T (f) is the unique mapping 
satisfying for all q G M 

x r(f){v(T,p )) = n(r,p ) = p(r,Po) - w(r,p ) = P v(t, po )(p( t iPo) ~ v ( r ,Po)) (5-49) 

where £>o = 9 + /(<?)• We can see X t {f) as the distance (in the sense of (|5.ip ) of the 
trajectory pt from M, starting at the time at a distance f from M. 

In the following, we will first prove that (|5.46p and (|5.47p imply that there exists an 
invariant manifold M £ for (|5.3p at distance e of M. Then we will prove (|5.46p and (|5.47p 
in three lemmas, paying attention on the relations between the different parameters. 

Suppose that the mappings X T satisfy (|5.46p for r G, [T, 2T] and that Xt satisfies (|5.47p . 
Then Xt has a unique fixed point in F(e, 1), which will be noted /o- Define <p £ = id + /o 
on M and M £ = <f>o(M). Since fo is a fixed point of Xt, if Po G Mq, then p^T G Mo for 
all k G N. Then to prove that M £ is an invariant manifold of (j5.3j) . it is sufficient to prove 
that for all t G (0, T), the functions /t defined by /t = X t (fo) are equal to /o- Using the 
property of semi-group and Xr(/o) = /o it is easy to see that ft = XT+t(fo), and thus 
(|5,46p implies that ft G J-(s, 1). But the same arguments show that ft is a fixed point of 
Xt for all t G (0,T). In conclusion, M e is invariant for (|5.3p . 

Now we prove ()5.46p and (|5.47p in the three following lemmas, which correspond to 
Lemmas 74.8, 74.9 and 74.10 in 



Lemma 5.4. There exists a £3 ^ E\ such that if e ^ £3, i/iere exists a 5^{e) of the form 
min(Ce, £3) such that if 5 53(e), we /ioue /or aZZ r G [T, 2T] and / G .F(e:, 1) 

||^r(/)||oo<e (5-50) 

Proo/ Let v G M, p = v + f(vo). We write (see (|5.49p ) 

XAf)(v(r)) = i* (r) (p(r) - vo) - P: {t) (v(t) - v ) ■ (5.51) 
The first term can be written as 

P: (t) ( P (t) - v ) = P: (T) (e-^o ( po _ VQ ) + F[p( S ) - Uo ] + G?bW]d«) • (5-52) 

Using the spectral gap, we bound the linear term 

\\P: {T) e- TLv °(p -v )\\ 2 ^ Cpse-^ T e (5.53) 

and the remaining term of (|5.52[) can be bounded in the same way as (|5.43p . Furthemore 
the second term of (|5.5ip is quadratic in e and 5, using a Taylor argument as in (|5.40p . 
Finaly, we get 

\\X T (f)(v T )\\ 2 ^ C A ({5 + e) 2 + 5))+Cpse-^ T e. (5.54) 
We supposed Cpse~ XlT ^ ^, thus we can choose 

e 3 = min(ei,j^-J and 5 3 {e) = muWei,£, — ej . (5.55) 



TRANSITIONS IN ACTIVE ROTATOR SYSTEMS 



25 



□ 

Lemma 5.5. There exists £4 ^ £3 such that if e ^ £4, there exists a 64(e) of the form 
min(Ce,e4) such that if 5 ^ 64(e), then for all f E J~(e, 1) we have X T (f) G T {e, 
for allr £[T,2T\. 

Proof It is sufficient to prove that X T (f) is Lipschitz with Lipschitz constant -gj- on all 
M n -B 2 (g, po) with po = 8C2e. Indeed in this case, if ||gi — q 2 \\ 2 > Po, then 

\\X T (f)( qi ) - X T (f)(q 2 )\\ 2 ^ -\\ qi - q 2 \\ 2 ^ -L.\\ qi - q 2 \\ 2 . (5.56) 

Po 4C 2 

Take ui,u 2 G M such that — u 2 \\ 2 ^ po an d / with Lipschitz constant I ^ 1. There 
exists viq,v 2 q G M such that Ui = v(r,pio) with p^o = ^io + f( v io)- Our goal is to show 
that under the hypothesis 

\\X T (f)( Ul ) - X T (f)(u 2 )\\ 2 \\X T (f)( Vl (r)) - X T (f)(v 2 (r))\\ 2 \\An(r)\\ 2 < 1 



IK -W2II2 hi(r) - v 2 (t)\\ 2 ||Au(t)|| 2 4C 2 ' 

(5.57) 

We use the decomposition 

An(r) = e- TL ^op^ o An(0) + An(r)-e- rL ^op^ o An(0) 

= e -^2 p* 2o An(0) + Ap(r) - Au(r) - e^^p^ Ap(0) + e^^p^ Av(Q) 
= e-^op; o A«(0) + Ap(r) - P° 20 Ap(i) + P° Ap(i) - A«(r) - e^of^ Ap(0) 
+ e- TL °™P° 20 Av(0) 

= [e-^o P S 2o An(0)] + [(P: 2(T) - P* 20 )Ap(r)] + [P* (Ap(r) - e^o A p(0))] 
+ [e-^ 2O p* 20 Av(0)] + [p° (r) Ap(r) - A,(r) 

We bound the first term using the spectral gap of L„ 20 , the second term using the smooth- 
ness of P s and Lemma [531 an d the third term in a similar way as (I5,43p , We use a Taylor 
decomposition for the two last terms, as in (I5.40p . Then we get (recall (I5.20P ) 

||An(r)|| 2 ^ (Cpse- XlT l + C 5 (p Q + 6 + e))\\Av(0)\\ 2 . (5.58) 



Since / is Lipschitz with Lipschitz constant I ^ 1 we have ||Ap(0)|| 2 ^ 2|| Au(0)|| 2 . Then 
using the part (3) of Lemma [531 we deduce 

||Au(0)|| a < 2||A«(t)|| 2 . (5.59) 

Furthemore we have chosen T such that Cpse~ XlT / 2 ^ yg, and thus Cpse~ XlT ^ 
(recall that C 2 = e Al/2 ). We obtain 

||Av(r)|| 2 8C 2 

Finally choose 

e 4 = min ( £ 3, 32 C 2 C 5 (l + 4C 2 ) ) = min ( £4 ' ,53 ( e )) ' ( 5 - 61 ) 

and the proof is complete. □ 
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Lemma 5.6. There exists £5 ^ £4 such that if £ ^ £5, there exists a 65(e) of the form 
mm(Ce, £5) such that for all f, L 6 T ( e, J : 



||*r(/i) - X r (/ 2 )|U < - / 2 |U . (5.62) 

Proof This time take i>io = u 2 o = v o and p^o = ^0 + /i(^o)- With the same decomposition 
as in Lemma 15.51 (with fewer terms, since v\q = V20) we get 

||An(T)|| 2 ^ (c P se- XlT + C 6 (<5 + e)) ||Ap(0)|| 2 . (5.63) 

We choose 

e 5 = min ( £ 4,Y^r^ and <^( e ) = min(e5,5 4 (e)) , (5.64) 
and in this case we get 

||An(T)|| 2 < -/ 2 [|oo- (5-65) 

Now notice that 

[|(Xt(/i)-*t(/ 2 ))MT))[|2 < HAnCrJIIz+llXTC/OCm^-XTC/OC^CT))^, (5.66) 
and since Xxifi) is Lipschitz with Lipschitz constant we get, using Lemma I5TT1 

||Xr(/i)(vi(T)) - X T (h)(v 2 {T))\\ 2 < ^||At;(r)|| 2 < - / 2 |U . (5.67) 

□ 

Proof of Theorem \2.1\ In these three lemmas, we see that if e is small enough, we can take 
5 proportional to e, thus adding a perturbation of type 5Cr[pt] to ()2. 1 j) creates an invariant 
manifold M$ situated at a distance 0(5) from M. It is proven in [24, (theorem 74.15, 
p. 531)] that the manifold M$ is C 1 in L\ and normally hyperbolic. Remark furthermore 
that (cf) S ) 1 (p) = v(p) for all p 6 M^. So to prove that <^> 5 is C 1 , it suffices to prove 
that v satisfies the hypothesis of the local inverse theorem between manifolds, that is Dv 
is a bijection between the tangent spaces of de two manifolds. Since the manifold is of 
dimension one, this property is implied by the lipschitz property of <p & . Furthermore we 
can estimate the differential of <j) & : ()5.60p for <p s gives an inequality for the local Lipschitz 
constant I s of 4>$ on all neighborhoods MLI B 2 (q, po) (po is introduced right before (|5.56p ): 

I 5 < ^r^ + CrS, (5.68) 

and we get that for a Cg > 

I 5 < C s 8, (5.69) 

which yields the bound we claim on the differential of <pg. □ 

Appendix A. On a norm equivalence 
The goal is to prove that the norms || • \\hi and || • \\y2 are equivalent, with 

- L q u := ^u" - [uJ * q + qj * u]' (A.l) 

and 

HI*? ■= \\(C + L q )u\\_ lil/q (A.2) 



TRANSITIONS IN ACTIVE ROTATOR SYSTEMS 



27 



with C > 0. Remark that by changing the constant C we get an equivalent norm. Since 



the norms 



-1,1/9 



and 



-i are equivalent, we will study \\(C + L q )u\\-y. We write 



inO 



u{6) = J2< 

J * q is of the type ae %e — ae~ ld , thus we can write 
Furthermore 



So if we denote 
then 

Consequently 
\\(C + L q )u\U 



2% 2% 

jo 



q{6) = 



qj * u 



Kay 
'^2A 



+1)0 



+ 



ifo- 



2/ 



1)5 



(A.3) 

(A.4) 

(A.5) 
(A.6) 
(A.7) 



6a n + n On - zan(a n _i - a n+1 ) + n— — c^_ x - n — - — c n+1 



(A.8) 



Suppose now that u G Hy. It is easy to see that there exists c > such that ||u||y2 ^ c||t/|| j h- 1 . 
Thus Yl n2 \ a n\ 2 < oo implies that ||(C + L 3 )u[|_i < oo and so Hy C V 2 . By expanding 
(TOl) and using Cauchy-Schwartz inequality we get 

|| (C + L,)u||_i > + (° 2 + n4 + 2Cn2 - «i n3 - a ^Cn) \a n \ 2 (A.9) 

where ay, a% ^ do not depend on u. It is clear that for C big enough ( depending on ay 
and 02 ) we have 



C 2 



1 



h n — avn > — n 

2 2 

r 2 

— + 2Cn 2 - a 2 Cn > 



(A.10) 
(All) 



and thus ||(C + 5? jIMI-Hi- We have shown that there exist c > such that for 

all u £ Hy, 

C^ 1 \\u\\y2 ^ ||u||if)_ ^ c||u||y2 . (A-12) 

But Hy is dense in V 2 (consider the finite sums of fourier series). If v £ V 2 , there exists a 
sequence v n in Hy such that v n -+ u for the VI 2 norm. Then v n is a Cauchy sequence for 
the Hy norm, and since Hy is complete, v £ Hy. In conclusion V 2 and Hy have the same 
elements. 

Remark A.l. By replacing (1 + n 2 ) -1 fry (1 + n 2 ) k , we can prove in the same way that 
|| (C + L q )u\\H k is equivalent to ||w||,ff fc+2 . Thus ||u||yn = || (1 + L ? ) n / 2 ti||_i ) i/ g is equivalent 

to IMItf n _i- 
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